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Background: Coronavirus envelope (CoV E) proteins have a predicted B-coil-B motif reported to target the Golgi complex. 
Results: This conserved domain forms B-structure on its own but is a-helical in the context of full-length SARS-CoV E protein. 


Conclusion: This domain is potentially involved in large conformational transitions. 
Significance: This is the first structural data of the extramembrane domain of any coronavirus E protein. 


Coronavirus envelope (CoV E) proteins are ~100-residue 
polypeptides with at least one channel-forming a-helical trans- 
membrane (TM) domain. The extramembrane C-terminal tail 
contains a completely conserved proline, at the center of a pre- 
dicted B-coil-B motif. This hydrophobic motif has been 
reported to constitute a Golgi-targeting signal or a second TM 
domain. However, no structural data for this or other 
extramembrane domains in CoV E proteins is available. Herein, 
we show that the E protein in the severe acute respiratory syn- 
drome virus has only one TM domain in micelles, whereas the 
predicted B-coil-B motif forms a short membrane-bound a-he- 
lix connected by a disordered loop to the TM domain. However, 
complementary results suggest that this motif is potentially 
poised for conformational change or in dynamic exchange with 
other conformations. 


Coronaviruses (CoV’; order Nidovirales, family Coronaviri- 
dae, subfamily Coronavirinae) are enveloped viruses organized 
into three groups (1, 2)*: group 1 (a-coronaviruses), group 2 
(B-coronaviruses) and group 3 (y-coronaviruses). Coronavi- 
ruses have been known to cause common cold symptoms in 
humans and a variety of lethal diseases in birds and mammals 
(4). However, in 2003, the virus responsible for the severe acute 
respiratory syndrome (SARS-CoV) (5) produced a near pan- 
demic with 8,273 cases and 775 deaths (6). In 2012, a novel 
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B-coronavirus (HCoV-EMC) (7-9) was discovered that has 
already led to many fatalities (10-12). 

The main coronavirus structural proteins are S (spike), E 
(envelope), M (membrane), and N (nucleocapsid), where §S, E, 
and M are integral membrane proteins. E proteins are ~100- 
residue-long polypeptides that are minor components in viri- 
ons but are abundantly expressed inside infected cells (13). 
They have a short hydrophilic N terminus, at least one pre- 
dicted terminal transmembrane (TM) domain, and a less 
hydrophobic C-terminal tail. 

Co-expression of E and M proteins is sufficient for formation 
and release of virus-like particles (14-19), and E proteins have 
been proposed to participate in inducing membrane curvature 
or in the scission of particles (20). Mutations at their C-terminal 
extramembrane domain impair viral assembly and maturation 
in the murine hepatitis virus (20). In the transmissible gastro- 
enteritis virus, the absence of E protein resulted in a blockade of 
virus trafficking in the secretory pathway and prevention of 
virus maturation (21, 22). In the case of the severe acute respi- 
ratory syndrome coronavirus (SARS-CoV), viruses lacking the 
E gene showed attenuation and did not grow in the central 
nervous system (23, 24), whereas stress response genes were 
up-regulated and cell apoptosis increased (25). These results 
suggested a role of E protein in both tissue tropism and patho- 
genicity, where modulation of stress responses contributes to 
viral attenuation. In fact, SARS-CoV AE attenuated viruses consti- 
tute promising vaccine candidates (26 —28). Thus, although E pro- 
teins are not absolutely essential for i vitro or in vivo coronavirus 
replication, their absence is clearly deleterious. 

Immunofluorescence experiments have shown that SARS- 
CoV E has acytoplasmically oriented C terminus and a lumenal 
N terminus (i.e. an N..oC cyto orientation) (13), which is consis- 
tent with the presence of a single TM domain. SARS-CoV E 
protein and other CoV E proteins have channel activity in syn- 
thetic membranes (29-31). This channel activity is mediated 
by formation of pentameric oligomers (32—34) and is only very 
mildly selective for cations (35). The only available structural 
data for CoV E proteins, obtained using synthetic TM peptides, 
is derived from the channel-forming TM domain in SARS-CoV 
E (32, 34, 36). No structural data are available for the predicted 
N- or C-terminal extramembrane domains despite the latter 
being critical for viral assembly (37, 38), although some results 
have been obtained using shorter synthetic peptides encom- 
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passing this domain (39). The C-terminal domain of E proteins 
contains a totally conserved proline residue, which in B- and 
y-coronaviruses is at the center of a predicted B-coil-B motif 
(Fig. 1A), reminiscent of viral internal fusion peptides (40 — 42). 
In SARS-CoV E, this motif was found to be responsible for 
redirecting a plasma membrane protein to the Golgi region. 
Conversely, mutations designed to increase its a-helical pro- 
pensity disrupted localization to membranes (43). 

The C-terminal tail of E proteins is also important for its 
interaction with the C-terminal domain of M protein (18, 44, 
45) at the cytoplasmic side of the endoplasmic reticulum-Golgi 
intermediate compartment, the budding compartment of the 
host cell. These interactions are the major drivers for envelope 
formation (46). The C-terminal tail of SARS-CoV E protein also 
interacts with the cellular protein PALS1 (47), which is thus 
depleted from the tight junctions in epithelial cells. Finally, 
SARS-CoV E interacts with the seven-domain (48) SARS non- 
structural protein 3 (Nsp3) (49). 

Thus, it is crucial to determine the structure of the C-termi- 
nal predicted extramembrane domain of E proteins, especially 
in the context ofa large construct that includes the TM domain. 
Until now, these structural studies have been hampered by 
problems in expression, purification, and stabilization of E pro- 
teins. We report herein the first detailed structure of a trun- 
cated form of the SARS-CoV E monomer that includes both its 
TM domain and its predicted B-coil-8 motif at the C-terminal 
tail. 


MATERIALS AND METHODS 


SARS CoV E Protein Constructs—Full-length SARS CoV E 
and a truncated version were used in this work. Full-length 
SARS-CoV E protein (E,;,) was cloned into pTBMalE plasmid 
downstream of the MBP fusion tag and tobacco etch virus 
cleavage sequence. The plasmid was transformed into Esche- 
richia coli strain BL21(DE3) codon plus for protein expression. 
The truncated form, Ey, was cloned into pNIC28-Bsa4 with an 
N-terminal His, tag followed by a tobacco etch virus cleavage 
sequence. The plasmid was transformed into E. coli strain 
BL21(DE3) Rosetta T1R for protein expression. In both con- 
structs, all three native cysteines (Cys-40, Cys-43, and Cys-44) 
were mutated into alanines. In addition, two E,,, mutants pre- 
viously described by Cohen et al. (43) were prepared by site- 
directed mutagenesis: (i) P54A (Eps4,) and (ii) V56A/Y57A/ 
V58A/Y59A (Egaza): 

Protein Expression and Purification—Non-labeled E protein 
was produced by growing the culture in Terrific broth medium 
at 37 °C until the culture density reached an Aggy of 2. Protein 
expression was induced by adding 0.5 mM isopropyl! 1-thio-B- 
p-galactopyranoside and growing the culture overnight at 
18°C. The cells were harvested by centrifugation at 7,500 X< g 
and stored at —80 °C. 

Stable isotope-labeled E protein was produced by growing 
the culture in LB medium at 37 °C. When the culture density 
reached an Ago, of 0.7, the medium was exchanged to M9 min- 
imal medium at 25% of the initial volume to achieve a high density 
culture, as described previously (50). The M9 medium was appro- 
priately supplemented with *"N-NH,Cl and [8C]glucose (Cam- 
bridge Isotope Laboratories) to produce '°N-labeled and 1°N/ 


12536 JOURNAL OF BIOLOGICAL CHEMISTRY 


C-labeled protein. Cultures were further grown for 1 h before 
inducing protein expression with 0.5 mM isopropyl 1-thio-B-p- 
galactopyranoside at 18°C. After 6 h, cells were harvested as 
described above and stored at —80°C. 

Frozen cell pellets were resuspended in lysis buffer (20 mm 
Tris, pH 8.0, 300 mM NaCl, 5 mn imidazole, 2 mm B-mercap- 
toethanol, and 10% glycerol) supplemented with 1 mm PMSF 
and 1.5% Triton X-100. The cells were completely lysed by son- 
ication and microfluidization. Insoluble particles were removed 
by centrifugation at 40,000 x g, and the supernatant was 
applied onto a pre-equilibrated nickel-nitrilotriacetic acid resin 
(Bio-Rad Profinity IMAC Ni?*-charged). The resin was washed 
with 20 mo Tris, pH 8.0, 300 mM NaCl, 20 mn imidazole, 2 mM 
B-mercaptoethanol, and 10% glycerol. Bound peptide was 
eluted in 20 mM Tris, pH 8.0, 300 mm NaCl, 250 mn imidazole, 
2 mM B-mercaptoethanol, 10% glycerol, and 5 mm myristyl sul- 
fobetaine (C14SB). E;p protein was directly TCA-precipitated 
and lyophilized, whereas E,,; was subjected to tobacco etch 
virus protease cleavage for 3 h at 30°C prior to precipitation 
and lyophilization. Further purification was achieved by using 
reversed-phase HPLC on a Phenomenex Jupiter C4 semi- 
preparative column (250 X 10 mm, 300-A pore size, 5-um par- 
ticle size). Lyophilized peptide was dissolved with 1% TFA in 
acetonitrile and separated under an isopropyl] alcohol-acetoni- 
trile linear gradient (4:1 (v/v) with 0.1% TFA). The identity and 
purity of peptide fractions were confirmed by SDS-PAGE and 
MALDI-TOF MS. 

Gel Electrophoresis—Standard SDS-PAGE was performed in 
13.5% Tris-glycine gel with TGS running buffer and stained 
with Coomassie Blue G-250. SDS-NuPAGE was performed in 
4-—12% NuPAGE® BisTris gel (Invitrogen) with NuPAGE® 
MES SDS running buffer and stained with SimplyBlue™ 
SafeStain (Invitrogen) according to the manufacturer’s protocol. 

To perform electrophoresis in the presence of perfluo- 
rooctanoic acid (PFO) detergent, we modified Invitrogen’s 
SDS-NuPAGE protocol by replacing SDS with PFO. Lyophi- 
lized peptide was dissolved in sample buffer containing 4% PFO 
and heated at 65 °C for 5 min prior to loading. The gel was run 
at 80 V for 2-3 h with MES running buffer containing 0.5% 
PFO. 

Blue native PAGE was performed as described previously 
(51). Lyophilized peptide was solubilized (0.1 mm) in sample 
buffer containing 25 mm SDS and either 25, 50, or 100 mm DPC. 
Aquaporin Z in 20 mmo SDS (heated at 65°C for 10 min) was 
included as an additional molecular weight marker. 

Fourier Transform Infrared Spectroscopy—Sample prepara- 
tion, data collection, and H/D exchange were performed essen- 
tially as described (32) on a Nicolet Nexus spectrometer (Mad- 
ison, WI). The peptides were incorporated in multilamellar 
liposomes by dissolving a dry mixture of 1,2-dimyristoyl-sn- 
glycero-3-phosphate (DMPC; Avanti Polar Lipids) and lyophi- 
lized peptide in HFIP at a 50:1 molar ratio. Fourier self-decon- 
volution was performed for some spectra using the following 
parameters: full width at half height, 20 cm~ '; narrowing factor, 
k = 1.5) (52). 

Analytical Ultracentrifugation—Sedimentation equilibrium 
experiments were performed using a Beckman XL-I analytical 
ultracentrifuge at 20 °C (53) and monitored by measuring the 
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absorbance at 280 nm. Lyophilized E;, peptides were dissolved 
at Aj gq of 0.3, 0.5, and 0.8 (12-mm path length cell) in 20 mm 
sodium phosphate, pH 5.5; 50 mm NaCl; and for detergent 
either 5 mM C14SB, 100 mm DPC, or 12.5, 25, or 50 mm SDS. To 
match the density of the SDS-DPC mixture, D,O was added at 
61.6, 65, and 72.4%, respectively, to each SDS concentration. 
The samples were centrifuged in six-channel charcoal-filled 
Epon centerpieces using quartz windows. A radial distribution 
profile was acquired after sufficient time to reach equilibrium, 
as tested by HeteroAnalysis. The data were processed and fitted 
to several monomer/n-mer models in SEDFIT and SEDPHAT 
(54). 

Circular Dichroism—CD data were acquired on a Chirascan 
CD spectrometer (Applied Photophysics) using a 0.2-mm 
quartz cuvette (Hellma). E+-p peptide samples were dissolved at 
1 mg/ml in 20 mM sodium phosphate, 50 mm NaCl, pH 5.5, 100 
mM DPC, with or without 50 mm SDS. CD spectra were 
acquired from 180 to 260 nm with a 1-nm spectral bandwidth 
and 3 replicates/spectra. Data points with excessive absorbance 
were excluded. After baseline subtraction, the data were ana- 
lyzed in Dichroweb (55) by using the CDSSTR method (56) and 
the SMP180 reference set (57). 

NMR Sample Preparation—Approximately 1.2 mg of lyoph- 
ilized E;-p protein was solubilized in 100 wl of methanol and 
dried under a dry stream of N, gas, resulting in a thin protein 
film deposit. The tube was placed in a vacuum lyophilizer over- 
night to remove any residual methanol. The thin protein film 
was then solubilized with sample buffer containing 20 mm 
sodium phosphate, pH 5.5, 50 mm NaCl, and 50 mm SDS. The 
sample was vortexed and sonicated several times until a clear 
solution was obtained, indicating protein reconstitution into 
detergent micelles. 

For paramagnetic relaxation enhancement (PRE) experi- 
ments, a single point mutation (S60C) was introduced into E+, 
by site-directed mutagenesis using appropriate sets of primers. 
Expression and purification protocol of the E;p-S60C mutant 
was the same as that of Erp protein. For labeling, 0.3 mm '°N- 
labeled E+-p-S60C was dissolved in 20 mm sodium phosphate, 
50 mM NaCl, 200 mm SDS, and 0.8 mm DTT at pH 5.5 and split 
into two equal portions for parallel labeling with (1-oxyl- 
2,2,5,5-tetramethyl-A*-pyrroline-3-methyl) methanethiosul- 
fonate (MTSSL) (TorontoResearch Chemicals Inc.) and a 
diamagnetic analog of MTSSL: (1-acetyl-2,2,5,5-tetramethyl- 
A®-pyrroline-3-methyl) methanethiosulfonate (AMTSSL; Toronto 
Research Chemicals Inc.). A 10-fold molar excess of both reagents 
was added from 75 mo stocks in methanol. The sample was 
vortexed for 30 min at high speed and incubated overnight at 
room temperature. A centrifugal filter unit (10,000 molecular 
weight cut-off; Millipore Corp.) was used to remove excess of 
both reagents. Labeled samples were washed four times by con- 
centrating to 100 pl. After a fourth wash, the sample was con- 
centrated to 180 pl for NMR measurements. 

Partial alignment of the E;., protein-micelle complexes rela- 
tive to magnetic field was obtained by using stretched poly- 
acrylamide hydrogels (58, 59). A 7% polyacrylamide gel was 
polymerized in a gel chamber of 5.4-mm inner diameter. After 
complete polymerization, gels were washed in H,O overnight 
and then twice with sample buffer containing 20 mm sodium 
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phosphate and 50 mm NaCl at pH 5.5. The gels were then com- 
pletely dried at room temperature. The protein solution con- 
taining E;,/SDS was soaked into the dried gels over 2 days to 
ensure complete rehydration. The hydrated 7% gel was then 
compressed into a 4.2-mm inner diameter open-ended tube 
using the gel press assembly (New Era Enterprise, Inc.). 

NMR Spectroscopy—NMR experiments were performed at 
308 K using an Avance-II 700 NMR spectrometer with cryo- 
genic probe. Sodium 2,2-dimethyl-2-silapentane-5-sulfonate 
was used as the internal reference for 'H nuclei. The chemical 
shifts of 7C and !°N nuclei were calculated from the 'H chem- 
ical shifts. The NMR data were processed using TopSpin 
version 3.1 and analyzed using CARA. Sequence-specific 
assignment of backbone 'HN, 1°N, '3C’, and '8C® was 
achieved by using two-dimensional H-'°N TROSY-HSQC, 
three-dimensional HNCO, HN(CA)CO, HNCA, HN(CO)CA, 
and HNCACB experiments on a '°N/13C-labeled E;, protein. 
Side-chain resonances were assigned using three-dimensional 
™N-resolved NOESY-HSQC (80-, 100-, and 150-ms mixing 
time), (H)CCH-TOCSY, and '°C-resolved NOESY-HSQC 
(120-ms mixing time). To identify membrane-embedded resi- 
dues, the NMR sample was lyophilized overnight and reconsti- 
tuted in 99% D,O. Immediately after reconstitution, two- 
dimensional 'H-'°N TROSY-HSQC was collected. For 
paramagnetic probe measurements, '7N-HSQC spectra were 
recorded in the presence and absence of 1 mM dry 5-doxyl ste- 
aric acid (5-DSA) or 16-doxyl stearic acid (16-DSA). Axially 
symmetric alignment tensor coefficients (axiality and rhombic- 
ity) were calculated using MODULE (60). The PRE effect was 
measured using '°N HSQC spectra of the S60C mutant before 
spin labeling and after MTSSL and dMTSSL labeling. The titra- 
tion experiments with HMA, Nsp3a, or SH(45—65) were per- 
formed with ‘°N-labeled E;z in a 1:4 molar ratio of SDS/DPC 
micelles at 318 K. Chemical shift perturbation (CSP) values 
were calculated using the following formula. 


CSP = ASH? + (0.23 x ABN)? (Eq. 1) 


Structure Calculation—NOE distance restraints were 
obtained from *N NOESY-HSQC (mixing time 80, 100, and 
150 ms) and 13C NOESY-HSQC (mixing time 120 ms) spectra, 
respectively. Backbone dihedral angle restraints (~ and ) were 
derived from 1°C’, °C%, 13C%, tH, and tH® chemical shift val- 
ues using TALOS+ (61). The short range and medium range 
NOE connectivities were used to establish the sequence-spe- 
cific 'H NMR assignment and to identify elements of the regu- 
lar secondary structure. Hydrogen bonds were derived from the 
H/D exchange experiment and NOE connectivity. 

Distance restraints were obtained from the measured PRE 
effect using the procedures described previously (62-64). The 
intensities of cross-peaks in the MTSSL (Jp) and dMTSSL (Id) 
were calculated in CARA. The correlation time was set to 10 ns. 
The ratios of intensities ([p/Id) were normalized against a set of 
the eight highest Jp/Id ratios, which were assumed to belong to 
peaks unaffected by PRE. For peaks with ratios below 0.15, no 
lower distance restraints were used, whereas upper restraints 
were set to 15 A. For peaks with ratios above 0.9, only upper 
restraints of 25 A were utilized. For peaks with ratios between 
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0.15 and 0.9, upper and lower distance restraints were gener- 
ated using +3 A margins. 

Structure calculations were performed using CYANA ver- 
sion 3.0 (65, 66) and visualized using PyMOL (Schroedinger 
LLC, New York). CNS 1.3 (67, 68) was used to refine the structure 
using the standard simulated annealing protocol. All of the 
restraints used in the calculations to obtain a total of 15 structures 
and all of the structure statistics are summarized in Table 1. 

Single Channel Activity Measurement—lIon channel activity 
of E,-, was measured by using Nanion Port-a-Patch®. Briefly, 
giant unilamellar vesicles of 1,2-diphytanoyl-sn-glycero-3- 
phosphocholine containing 10% cholesterol were prepared in 
1 M sorbitol using Nanion Vesicle Prep Pro®, following the 
manufacturer’s protocol. Into 100 yl of a giant unilamellar ves- 
icle solution, 0.5 pl of Ey, peptide solution in ethanol (0.1 
mg/ml) was added and incubated for 1 h at room temperature. 
E;r-containing giant unilamellar vesicles were subsequently 
deposited onto 6—8 megaohm NPC®-1 chips (Nanion). Con- 
ductance was measured under symmetrical buffer conditions 
(10 mm HEPES, 500 mm NaCl, pH 5.5). 

Surface Plasmon Resonance—The nsp3a sequence was sub- 
cloned from pcDNA3(+) into pET28b upstream of a C-termi- 
nal His, tag for expression in E. coli. The protein was expressed 
and purified as described previously by Serrano et al. (69). A 
negative control, consisting of C-terminal peptide from the 
small hydrophobic (SH) protein of human respiratory syncytial 
virus (RSV SH(45-—65)), was synthesized by standard solid 
phase and purified by reverse-phase HPLC. Surface plasmon 
resonance measurements were performed on a Biacore 3000 
system (GE Healthcare) using 10 mm phosphate buffer at pH 
6.5, 100 mM NaCl, 3mm EDTA, 0.05% n-octyl- B-p-glucopyran- 
oside, and 0.27% C14-betaine at 25 °C. E,., was immobilized to 
15,000 RU onto a research grade CM5 sensor chip (GE Health- 
care) using standard amine-coupling chemistry. Briefly, a buf- 
fer-equilibrated carboxymethyl dextran surface was activated 
with a 10-min injection of a 1:1 mixture of 0.05 m N-hydroxy- 
succinimide and 0.2 m N-ethyl-N-[3-(diethylamino)propyl]car- 
bodiimide. E--, peptides dissolved in 10 mM sodium acetate, 15 
mM DPC (pH 5.0) were passed over the activated surface to 
achieve the desired response level. Another 10-min injection of 
1 Methanolamine-HCl (pH 8.5) was used to deactivate the sur- 
face and remove any non-covalently bound protein. Kinetic 
measurements of immobilized E+, association with Nsp3a and 
SH(45— 65) (49 nm to 25 uM in 10 2-fold serial dilutions) were 
performed with a 1-min association phase and 5-min dissocia- 
tion phase at a 30 pl/min flow rate. Each concentration was 
tested in duplicate. No regeneration was necessary because all 
complexes dissociated within the monitored time. Sensor- 
grams were double-referenced (70) and globally fit to a steady- 
state model to obtain affinity values. 


RESULTS AND DISCUSSION 


Expression and Purification of SARS-CoV E Protein—Ini- 
tially, we successfully expressed and purified full-length SARS- 
CoV E protein (E,,, Fig. 1B) by using either a B-barrel (71) or 
MBP as fusion tags. However, the yield of pure protein was low 
due to the presence of truncations (not shown). Nevertheless, 
the pure sample obtained (Fig. 1, C and D) was sufficient for 
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backbone assignment in SDS micelles, although not to produce 
a three-dimensional model. Therefore, a series of hexahisti- 
dine-tagged SARS-CoV E constructs were screened to obtain 
an expressing and well behaved sample. The best construct, 
encompassing residues 8—65, was successfully expressed in 
E. coli and purified by affinity chromatography in milligram 
amounts without any enzymatic cleavage steps. This truncated 
construct (Ep; see Fig. 1B) has an N-terminal His tag and a 
16-residue linker that connects it to residues 8-65. The puri- 
fied peptide appeared in MALDI-TOF MS as a single-charged 
peak at 8,997 Da and a double-charged peak at 4,512 Da, con- 
sistent with the calculated molecular mass of Eyp, 8,995 Da, and 
a small proportion of larger (dimer to pentamer) oligomers (Fig. 
1E). After HPLC, the E;-g monomer (9 kDa) showed anomalous 
migration in standard SDS-PAGE (Fig. 1F), as reported previ- 
ously for the full-length protein (31). 

Identification of the Membrane-embedded Region of Erp by 
NMR—Screening of reconstitution conditions identified SDS 
as the best environment to achieve good peak dispersion in 
both dimensions (Fig. 2A). Comparison of 'H-'°N TROSY- 
HSQC spectra in water and in 99% D,O (Fig. 2, A and B) iden- 
tified a stretch of 21 residues, from Leu-18 to Leu-39, protected 
from H/D exchange. The same residues were protected from 
H/D exchange in E,, (not shown), which indicates that both 
Exp and E;, have only one TM domain. 

The topology of E,, was also delineated by introducing the 
paramagnetic probes 5-DSA and 16-DSA. Upon the addition of 
5-DSA, a slight intensity reduction was observed around resi- 
dues 11-20 and 40-55 (Fig. 2C), suggesting that these residues 
are located at or near the micelle surface. The addition of 
16-DSA caused pronounced intensity reduction in the stretch 
of TM residues 19-40 and periodically in the stretch 48-61 
(Fig. 2D). Together with the H/D exchange data (Fig. 2, A and 
B), these results again demonstrate the presence of a single TM 
domain and suggest that the stretch of residues 55-65 corre- 
sponds to a domain bound to the surface of the micelle. Finally, 
a ‘H-1°N steady-state heteronuclear NOE experiment (Fig. 2E) 
showed that the protein forms a well folded structure, with high 
™H/°N NOE values. However, more dynamic regions (lower 
™H/1°N NOE values) can be observed (i) at the N terminus and 
(ii) at residues 46—55, which connects the TM domain to the 
proposed membrane-bound domain. The latter is predicted to 
be a-helical (Fig. 2F). 

NMR Structure Determination of Ep in SDS Micelles; Struc- 
ture Calculation—The restraints shown in Table 1 were used in 
a calculation to obtain a total of 15 structures, with a root mean 
square deviation of 0.27 + 0.11 and 0.70 + 0.13 A for backbone 
and all heavy atoms, respectively (Fig. 3A). A longer a-helix 
(residues 15-45) encompasses the TM domain, which is con- 
nected to another shorter C-terminal a-helix (residues 55-65) 
by a flexible (see Fig. 2£) linker domain (residues 46 —54), form- 
ing an L-shape. The short extramembrane helix may be partially 
bound to the micelle surface, as suggested by the pattern of inten- 
sity attenuation of paramagnetic reagents (see Fig. 2, B and C). The 
most affected residues, which would face the micelle surface, are 
Val-52, Thr-55, Tyr-59, and Lys-63 (Fig. 3B). 

The TM a-helix has a slight bend at residues 26 —30, consis- 
tent with previous results obtained for the synthetic TM 
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FIGURE 1. Sequences, expression, and purification of SARS-CoV E,,. A, alignment of representative sequences of E proteins in a-, B-, and y-coronaviruses. 
The cysteine residues are underlined, the conserved proline is highlighted (gray), and the four residues mutated to alanine in the E,,,, mutant (see “Materials 
and Methods”) are shown in red. For these four proteins, the prediction of secondary structure is shown below in a color code, with the TM domain indicated as 
a black line; B, proteins used in the present work: a His-tagged construct (E;z) encompassing residues 8 - 65 (boldface type, underlined), and full-length SARS-CoV 
E (E,,). In E-,, the fragment SNA results from the cleavage of the tag. In both proteins, the native cysteines were mutated to alanine (C40A, C43A, and C44A; see 
asterisks); Cand D, MALDI-TOF MS spectra (C) and standard SDS-PAGE (D) of pure E,, with the species labeled; E and F, same as for purified E;,; the identities of 
various single- and double-charged species are indicated. The calculated mass of E;p is 8,995 Da. 
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TABLE 1 
Restraints and structure statistics for the selected 15 structures of E,, 
NMR restraints 
Total unambiguous distance restraints 2,258 
Intraresidual 1,334 
Sequential (|i — j] = 1) 230 
Short-range (|i — | < 1) 1,564 
Medium (2 = |i — j| = 4) 321 
Long range (|i — j| = 5) 11 
Dihedral angle restraints 85 
Hydrogen bond restraints* 20 
RDC restraints 44 
PRE restraints 38 
Root mean square deviation from the experimental 
residual dipolar couplings (Hz) 
1D ee 0.71 + 0.03 
Root mean square deviation from the average 
atomic coordinates (residues 12-63, A)’ 
Backbone atoms 0.27 + 0.11 
All heavy atoms 0.70 + 0.13 
Ramachandran analysis (%) 
Residues in most favored regions 87.5 
Residues in additional allowed regions 12.5 
Residues in generously allowed regions 0.0 
Residues in disallowed regions 0.0 


# Backbone hydrogen bonds of a-helix were applied to regions confirmed to be 
a-helical, according to the local NOE pattern and HN-H,O chemical exchange 
experiments. 

° Statistics were calculated and averaged over an ensemble of 15 structures with 
lowest target function according to CYANA. 


domain in DPC micelles (34). A kink near this location is sug- 
gested by the short distance (2.0 + 0.1 A) between Thr-30 Hy 
and the carbonyl oxygen at Phe-26, in the range of a hydrogen 
bond (Fig. 3C). 

Effect of Truncation and Environment on E;p Secondary 
Structure—To assess the effect of the truncation and the pres- 
ence of a His tag on the E;, secondary structure, E,, and E,, 
were compared. The possible effect of the reconstitution envi- 
ronment was also determined. 

CD and IR Spectra of Erp and E,, in Detergent and Lipid 
Membranes—The CD spectra of E-p in DPC, SDS, and mixed 
(1:2 molar ratio) SDS/DPC micelles are almost superimposable 
(Fig. 4A) with minima at 209 and 222 nm. Also, these data are 
entirely consistent with the CD spectra of E,,, which was pre- 
dominantly a-helical in both SDS and DPC micelles (71). When 
reconstituted in DMPC membranes, Ey, and E;,; produced an 
almost identical spectrum, with an amide I band centered at 
1,655 cm ' (Fig. 4B) characteristic of a predominantly 
a-helical conformation. Overall, these data show that both 
Exp and Ep, (i) are predominantly a-helical and (ii) have a 
secondary structure that is not significantly affected by the 
reconstitution environment, supporting the relevance of the 
Exp structure (Fig. 3). 

?3Ca Chemical Shifts—'*Ca chemical shifts are highly corre- 
lated with secondary structure (72, 73). Comparison of *Ca 
chemical shifts of Ep and E;; in SDS (Fig. 4C) shows that E;p 
shifts (red) are almost identical to their counterpart residues 
(residues 8—65) in E,, (blue). In E,,, the shifts for the last 10 
residues (residues 66 —76) are predicted to correspond to ran- 
dom coil, which is consistent with results obtained previously 


Structure of a Coronavirus Envelope Protein 


with a synthetic peptide spanning residues 59-76 (39). This 
peptide produced a broad amide I band in the IR spectrum 
centered around 1,650 cm ~' and immediately experienced 
complete H/D exchange. 

When the *°C shifts for Eg were compared for SDS and an 
SDS/DPC mixture (1:4 molar ratio), only the N-terminal region 
(residues 7-15) showed minor differences. Finally, a compari- 
son between E,, and E,;z in an SDS/DPC mixture (1:4 molar 
ratio) also revealed an almost identical pattern in the C-termi- 
nal tail (not shown), although the spectral resolution for E,, 
was reduced in other regions. Overall, these results are consis- 
tent with those shown above (Fig. 4, A and B), indicating that 
the detergent used, the tag, and the truncation has a minimal on 
Ep secondary structure. 

Ion Channel Activity of E7p and Its Inhibition by HMA—Pu- 
rified E,,; has channel activity (71), which is inhibited by the 
drug HMA (30). An J/V plot obtained in a symmetrical 0.5 Mm 
NaCl experiment for purified E,, in 1,2-diphytanoyl-sn- 
glycero-3-phosphocholine (Fig. 54) was used to determine a 
conductance of 0.39 + 0.02 nanosiemens. For comparison, syn- 
thetic full-length SARS-CoV E and E,,, (residues 7—38) pro- 
duced single channel conductances of 0.19 + 0.06 picosiemens 
and 0.18 + 0.12 nanosiemens in 1 M NaCl (35), although the 
values were higher in 1 M KCI, with 0.37 + 0.16 and 0.31 + 0.12 
for full-length SARS-CoV E and E,,,,, respectively. The lower 
conductance observed in synthetic samples may be due to 
extraneous modifications or impurities resulting from expo- 
sure to harsh chemicals. Representative traces of E;, channel 
activity (Fig. 5B) and complete inhibition after the addition of 
10 um HMA (Fig. 5C) suggest that Ep is entirely functional. 

To determine the binding site for HMA, we measured the 
differences in CSP values before and after the addition of the 
drug (Fig. 5, D and £). In SDS, the average CSP value was low, 
0.006 ppm (Fig. 5F), even at an HMA/E+p 10:1 molar ratio, 
suggesting no significant binding. However, the same panel 
shows that the addition of HMA to E.,, in mixed SDS/DPC (1:4 
molar ratio) micelles produced an average CSP value of 0.013 
ppm, even at an HMA/E-p 2:1 molar ratio. 

Except for Asn-64 and Leu-65, the residues that showed sig- 
nificant CSP (CSP = 0.025 ppm) clustered near the membrane 
interface regions of the TM domain (Fig. 5F, see arrows). Near 
the N-terminal side of the TM domain, the most affected were 
Glu-8, Gly-10, Thr-11, Val-14, Asn-15, and Ser-16, the latter 
two consistent with observations made on the TM channel in 
DPC (34). At the C-terminal end of the TM, Leu-37 was the 
most affected, suggesting that the interaction of HMA at He of 
Arg-38 reported previously (34) may have been an artifact due 
to the use of a TM peptide. 

Oligomeric State of Err 

Gel Electrophoresis—The localized changes in chemical 
shifts observed after HMA addition to E;p in mixed DPC/SDS 
micelles (Fig. 5£), but not in SDS micelles (Fig. 5D), suggest that 
binding sites for HMA may have been induced after DPC addi- 


FIGURE 2. Topology and secondary structure of E,g. 'H-'°N TROSY-HSQC spectra of 0.2 mm E;p in 50 mm SDS in H,0 (A) and in 99% D,0 (8). The cross-peaks 
are labeled by one-letter code and residue number; C and D, peak intensity reduction upon the addition of 5-DSA (C) and 16-DSA (D), calculated as the ratio of 
peak intensity before and after the addition of the paramagnetic reagents; £, 'H-'°N steady-state heteronuclear NOE experiment; F, sequential and medium- 
ranged NOE connectivity between residues, displayed as bands under the respective residues. 
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FIGURE 3. Structural model of E,,. A, superposition of an ensemble of 15 calculated simulated annealing structures of E,, (only the sequence corresponding 
to E protein, 8-65, is shown). Side chains are shown as line representations; the residues at the ends of the two helical segments are indicated. B, residues of the 
C-terminal extramembrane a-helix oriented toward the micelle surface (blue). C, ribbon representation of the TM central region, with the carbonyl oxygen of 


Phe-26 forming a hydrogen bond to the side chain of Thr-30. 
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FIGURE 4. Equivalence in secondary structure of E,, and E,,. A, CD spectra of E;, in DPC (black), 1:2 molar ratio SDS/DPC mixture (blue), and SDS (red). B, 
infrared amide | band of E;p (red) and E,, (b/ue) in DMPC lipid bilayers and their respective Fourier self-deconvolved spectra (dotted lines). C, comparison of 
secondary '2Ca chemical shifts (deviation from tabulated random coil '*Ca chemical shift values) for E;p (red dots) and E,, (blue dots) in SDS micelles and for Ex, 
in (1:4 molar ratio) mixed SDS/DPC micelles (white dots). For the latter, Pro-54 and Thr-55 were excluded from the analysis due to significant line broadening; 


Arg-38 was excluded from the analysis due to the peak overlapping. 


tion to SDS, possibly through E;, increased oligomerization 
and population of the pentameric form. Thus, we examined the 
oligomeric state of E,-, using gel electrophoresis performed in 
the presence of SDS or an SDS/DPC mixture. In contrast with 
the results obtained in normal SDS-PAGE (Fig. 1), we used 
NuPAGE, where E,, migrated with its expected monomer 
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molecular weight, as a single band at 9 kDa (Fig. 6A). No oligo- 
merization is consistent with the lack of interaction between 
HMA and E,, observed in SDS micelles. 

To test the effect of DPC, in a blue native polyacrylamide gel, 
a constant SDS concentration (25 mm) was titrated with 
increasing DPC (Fig. 6B), from 1:1 to a 1:4 molar ratio. Bands 
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FIGURE 5. Channel activity of E,, and interaction with HMA. A, //V plot for E;, in 1,2-diphytanoyl-sn-glycero-3-phosphocholine bilayers in a symmetrical 
experiment where both cis and trans compartments contained 10 mm HEPES and 500 mm NaCl at pH 5.5. Each point represents the mean of at least three 
current readings. The /ine is a linear regression fit of data points, which produced a slope of 0.39 + 0.02 nanosiemens. B, selected traces of 12 s each, recorded 
at various holding potentials of E;g. C, channel activity recorded at 60 mV holding potential and after the addition of 10 zm HMA (arrow). D, 'H-'°N TROSY-HSQC 
spectra of 0.1 mm '°N-labeled E,g in SDS micelles (D) before (b/ue) and after (red) the addition of 1 mm HMA. E, same as for 0.2 mm '°N-labeled E;p in (1:4 molar 
ratio) SDS/DPC micelles upon titration with 0.4 mm HMA. Some shifts are indicated with arrows; F, CSP of the backbone amide resonances of E;, before and after 
the addition of HMA in SDS (red) and (1:4 molar ratio) SDS/DPC micelles (blue). Note that the HMA/E;g molar ratio was 10 in SDS and only 2 in SDS/DPC micelles. 
The arrows show residues with significant change in chemical shifts after the addition of HMA. The TM domain is indicated only to guide the eye. 


consistent with dimers and trimers were observed at 1:1 and 1:2 
molar ratios, whereas tetramers and pentamers were observed 
at a 1:4 molar ratio. This increasing oligomerization is again 
consistent with the binding of HMA observed in the presence of 
DPC. 

Sedimentation Equilibrium, SDS/DPC Micelles—For a sam- 
ple corresponding to a 1:4 molar ratio of SDS/DPC, the equilib- 
rium sedimentation data (Fig. 7A) could be fitted to several 
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oligomeric models, from trimers to heptamers (Fig. 7B). Similar 
ambiguous results were obtained for a 1:2 molar ratio SDS/DPC 
mixture (Fig. 7C), suggesting that E,., forms a mixture of 
oligomers in SDS/DPC micelles, consistent with electrophore- 
sis results (Fig. 6). The fit of the traces (Fig. 7A) to a monomer- 
pentamer model produced an apparent K, of 1.25 X 101° m + 
(i.e.6 X 10°m +) ora K, for monomer-monomer interaction of 
0.17 mM and a molar fraction standard free energy change (74), 
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FIGURE 6. Gel electrophoresis of E,, in SDS and in a SDS/DPC mixture. 
Shown are Exp in 44-12% SDS-NuPAGE (A) and 4-16% blue native PAGE (B) of 
Erp in 25 mm SDS with increasing concentration of DPC, as indicated. E. coli 
aquaporin Z (AQPZ) was included as an additional molecular weight marker. 
Bands and oligomeric states are indicated by arrows and black dots, 
respectively. 


A 10, 


AQPZ 


Ase 


AG,°, of approximately —4 kcal‘mol~*. Overall, the above 
results show that higher oligomeric states, including pentam- 
ers, are observed when increasing the DPC concentration, con- 
sistent with the larger shifts observed after adding HMA to 
mixed DPC/SDS micelles. 

Ey Forms Pentamers in C14 Betaine and PFO—We have 
shown previously that SARS-CoV E TM domain and E,,;, form 
pentamers (36, 71). When E.-p was solubilized in C14 betaine 
detergent (Fig. 7D), data could be optimally fitted to a mono- 
mer-pentamer model (K, = 10'°m *) witha significantly lower 
x’ than that obtained in SDS/DPC. Also, the mobility of Ez, in 
PFO-NuPAGE (Fig. 7E) corresponds to a molecular mass of 45 
kDa (i.e. a pentameric form). Thus, E-TM, E+, and E,; show a 
similar oligomerization behavior, being able to form pentam- 
eric channels that are inhibited by HMA. 

Interaction of Ep with Nsp3a—The interaction of Nsp3 with 
SARS-CoV E protein, thought to be related to E ubiquitination, 
was previously mapped to an N-terminal acidic domain, Nsp3a 
(49). Although the site of interaction with E is not known, it is 
likely to involve the C-terminal tail because it is the largest 
extramembrane domain. Thus, to test that the conformation of 
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FIGURE 7. Sedimentation equilibrium of E;, in SDS/DPC micelles. A, radial distribution profiles (open circles) of Erg in a (25:100 mm) SDS/DPC mixture at 
16,600 rpm (red), 20,300 rpm (green), and 24,900 rpm (blue). The profile was fitted to a monomer-pentamer self-association model (black line), and the fitting 
residuals are shown below each plot. B-D, global reduced y? values obtained after data fits to different monomer/n-mer models of E;, association in SDS/DPC 
micelles (B), a 50:100 mm SDS/DPC mixture (C), and 5 mm C14-betaine (D). EF, 4-12% PFO-NuPAGE of E;p. 
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FIGURE 8. Interaction between E,, and Nsp3a. A, sensorgrams corresponding to the interaction between purified RSV SH(45-65) and immobilized E;p (red). 
The steady-state model (dark red) did not produce a good fit toa 1:1 model of interaction. The association phase extends from 0 to 60 s, whereas the dissociation 
phase extends for minutes. B, same as for Nsp3a and immobilized E;, (blue) and fit to a steady-state model (red). Inset, dose-response plot, where the 
equilibrium responses of Nsp3a were plotted against the log,,(concentration) of Nsp3a. Although the fit (b/ue line) yields an affinity of 1.6 mm, binding is already 
evident even in the interval 1-10 wm Nsp3a concentration. C, CSP of the backbone amide resonances of 0.2 mm '°N-labeled Ey, in (1:4 molar ratio) SDS/DPC 
micelles upon titration with 0.4 mm Nsp3a (blue dots) or the negative control RSV SH(45-65) (red dots). 


Exp is structurally equivalent to the corresponding sequence in 
E;1, we tested its ability to interact with Nsp3a (49) using sur- 
face plasmon resonance. 

As a negative control, the mock peptide RSV SH(45- 65) (see 
“Materials and Methods”) was used to bind to immobilized E;, 
in detergent (see “Materials and Methods”) (Fig. 8A). Binding 
and unbinding of RSV SH(45-—65) was very slow (>1 min), and 
the data could not be fitted to any model, suggesting a nonspe- 
cific interaction. 

In contrast, rapid and reversible interaction was observed 
using purified Nsp3a (Fig. 8B). These data could be fitted with a 
stoichiometry of 1:1, although the end point could only be esti- 
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mated. From that model, affinity was determined as 1.6 mm, 
although we note that binding is already evident at concentra- 
tions as small as 1-10 po (see Fig. 8B, inset). 

The residues involved in the interaction between Nsp3a and 
Ep were identified by the differences in CSP values before and 
after the addition of Nsp3a to E;, in SDS/DPC micelles (1:4 
molar ratio). Large chemical shift changes at Leu-39, Val-49, 
and Leu-65 (Fig. 8C) indicated a potential binding site at the 
C-terminal tail, whereas the central TM region did not experi- 
ence any change. Overall, these data suggest that E,, has a 
structure in SDS/DPC micelles that is similar to the native fold 
of E,, in biological membranes. 
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FIGURE 9. Structural flexibility at the putative B-coil-6 motif. A, amide | band corresponding to E,,, Eps4q, aNd E,,, 4 in DMPC bilayers. Regions that change 
after mutation of the four residues indicated in Fig. 1A (E4,, 4) are shown as arrows. B, Fourier self-deconvolved spectra corresponding to the amide | bands 
shown in A. C, infrared amide | and II bands corresponding to SARS-CoV E peptide E4¢_ 6 in DMPC lipid bilayers before (b/ue) and after (red) being exposed to 
D,0 (39). D, possible equilibrium between two conformations at the putative B-coil-B motif, shifted toward an a-helical form, between the model determined 
experimentally for E;z and E,, and one built with prediction tools (PEP FOLD) (3) and consistent with data obtained with synthetic fragment E4¢_ 6. shown in C. 


FIGURE 10. Pentameric model formed by E,,. Side (A) and top (B) views of a proposed E;, pentamer structure shown in ribbon representations. The side chains 
of Val-49 and Leu-65, which have been shown to interact with HMA, are shown as line representations. 


Structural Discordance at the C-terminal Tail—Having 
determined the suitability of the construct used by comparison 
with the full-length protein and shown that the nature of the 
detergent did not affect the secondary structure significantly, 
we investigated the structure of the C-terminal extramembrane 
domain of E,, and its apparent contrast with secondary struc- 
ture predictions. 

The secondary structure predictions clearly suggest the pres- 
ence of a B-turn-B motif in the C-terminal region, not only in 
SARS-CoV E but also in other representative sequences of 
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other coronaviruses (see Fig. 1A). Indeed, we have reported 
previously that the synthetic peptide E,, <9, which encom- 
passes this putative B-hairpin, folds as B-strands, is completely 
resistant to H/D exchange, and has a very high tendency to 
aggregate in solution (39). However, we have shown here that in 
the context of E,, or Ep, this domain does not adopt a B-struc- 
ture. To test if at least some part of the population adopts this 
structure in lipid membranes, we mutated E,, residues 56—59 
to increase the predicted helicity of this part of the molecule 
(43), from B-branched Val or bulky Tyr to small side chain Ala, 
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to obtain the construct E,,;,. Last, to test the effect of the 
conserved Pro-54 at the center of this putative motif, we gen- 
erated the mutant Eps,,. 

The IR spectrum of these mutants when reconstituted in 
DMPC membranes (Fig. 9A) shows a similar amide I band. 
However, for the mutant E,,; 4, two shoulders are eliminated 
(see arrows). This indicates that these shoulders in E,, do not 
represent misfolded protein, but they may correspond to the 
bona fide B-structure conformation present in a small part of 
the population. However, the significance of these mutations is 
not completely clear because, although the mutations intro- 
duced in E44; prevented Golgi complex accumulation when 
the C terminus tail of SARS-E was coupled to VSV-G (43), a 
similar effect was also observed for the P54A mutant, which in 
our IR spectrum showed no obvious differences with respect to 
Exp 

As stated above, the synthetic peptide E,, ¢, (39), which 
encompasses the predicted B-hairpin in SARS-CoV E, was 
found to produce ~100% f-structure and was completely 
resistant to H/D exchange (Fig. 9C). Combined with the effect 
observed for the E,,,; , mutant, we propose that this B structure 
may be in dynamic equilibrium with the much more abundant 
a-helical form (Fig. 9D). A delicate balance between these two 
forms may alter processes in the infected cell (e.g. membrane 
scission, binding to protein partners, or E protein localization). 

Finally, the HMA titration results showed large CSP values 
for residues Val-49 and Leu-65, which are far apart in the 
sequence. Using the SymmDock server (75, 76), a reconstructed 
E;p pentameric model was obtained based on the published 
structure of the E-TM pentamer (34) and the current E.-, struc- 
ture. The model suggests that these two residues may be in fact 
spatially close (Fig. 10) and belong to different monomers, pro- 
viding a rationale for the observation above. 

CoV E proteins have been proposed to have at least two roles. 
One is related to their TM channel domain. This would be 
active in the secretory pathway, altering lumenal environments 
and rearranging secretory organelles and leading to efficient 
trafficking of virions (38, 77). The other would be related to 
their extramembrane domains, particularly the C-terminal 
domain. This is involved in protein-protein interactions and 
targeting, among other roles. 

E protein participates in M-M and E-M interactions (17, 44), 
which are interesting targets for drug discovery. Also, forma- 
tion of viral particles appears to be facilitated by a broad range 
of E sequences (78), which suggests that a common topology is 
more important than sequence requirements. E proteins have 
been suggested to act as chaperones during packaging (79), but 
the precise mechanism by which this takes place is not known. 
In this context, the structure determined here sheds light on a 
critical domain present in most CoV E proteins. 
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